Electromagnetically Induced Transparency in strongly interacting Rydberg Gases 
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We develop an efficient Monte-Carlo approach to describe the optical response of cold three- 
level atoms in the presence of EIT and strong atomic interactions. In particular, we consider a 
" Rydberg- EIT medium" where one involved level is subject to large shifts due to strong van der 
Waals interactions with surrounding Rydberg atoms. We find excellent agreement with much more 
involved quantum calculations and demonstrate its applicability over a wide range of densities and 
interaction strengths. The calculations show that the nonlinear absorption due to Rydberg- Rydberg 
atom interactions exhibits universal behavior. 



The effect of electromagnetically induced transparency 
(EIT) jl[ in light-driven multi-level systems continues to 
play a pivotal role in quantum and nonlinear optics. En- 
abling slow light propagation and thus long photon inter- 
action times at low loss levels 0, Q , EIT media provide a 
promising route to applications in optical communication 
and quantum information science. Optical nonlinearities, 
however, typically arise from higher order light-atom in- 
teractions, such that realizations of such applications at 
very low light intensities |3| remain challenging. In this 
respect, recent experimental studies of EIT in cold Ry- 
dberg gases [E-Q are opening up new perspectives for 
nonlinear optics on a few photon level. Exploiting the 
exaggerated properties of Rydberg atoms and, in partic- 
ular, the strong interactions among the atoms, nonlinear 
phenomena can be greatly enhanced in ultracold Ryd- 
berg gases. 

Recently, several different methods have been used to 
study laser-driven interacting gases (lol - [T^ . A simul- 
taneous treatment of EIT and long-range interactions, 
however, poses additional challenges. Common meanfield 
approaches, as successfully applied to radiation trapping 
effects on EIT are found to fail 11] due to the large 
strength and enormous range of the van der Waals in- 
teraction, which lead to non-negligible correlations be- 
tween the atoms. On the other hand, exact descriptions, 
based on Hilbert space truncation by interaction-blocked 
many-body states [H, [11], are also inapplicable since 
there always remains an exponentially large number of 
many-body states involved in the interaction-free probe 
transition. 

Here we present a theoretical approach that allows to 
obtain the fully correlated steady state populations via 
classical Monte Carlo sampling. This is shown to yield 
the nonlinear optical response to classical light fields in 
the presence of arbitrarily strong atomic interactions. A 
comparison to reduced density matrix calculations shows 
very good agreement for small and moderate densities. 
Upon proper scaling the simulation results reveal a uni- 
versal behavior of the nonlinear absorption, which illus- 
trates the role of Rydberg atom interactions in the emer- 
gence of dissipative photon-photon interactions. 

The considered Rydberg-EIT level scheme is shown in 
FiglH The signal laser couples the ground state |1) to a 



FIG. 1. Illustration of the considered three-level ladder 
scheme. Two laser fields successively couple the states |1), |2) 
and |3). For isolated atoms, EIT is realized on two-photon 
resonance Ai — — A2. The van der Waals interaction shift 
Ce/rfj between atoms in the Rydberg state |3) modifies this 
ideal transparency and leads to a nonlinear optical response 
of the medium. 



low lying state |2) with Rabi frequency fli. State |2) is 
coupled by a strong control laser {^2 > fii) to a highly 
excited Rydberg state |3). The resulting dynamics of 
a gas composed of N such independent atoms is, thus, 
governed by the Hamiltonian 
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Ho = --^(Ai|2,)(2,| -I- (Ai + A2)|3,)(3,| 
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In addition, the intermediate state |2) radiatively de- 
cays with a rate 7, whereas spontaneous decay of the 
long-lived Rydberg level can safely be neglected. For 
resonant driving, each atom settles into a dark state, 
\di) ~ ri2|li) — ^^i|3i), which is immune to the laser cou- 
pling and radiative decay [ij . Consequently, the complex 
susceptibility xi2 of the lower transition vanishes and the 
medium becomes transparent. 
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In the presence of van der Waals interactions 
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between atoms in state |3) the situation becomes more 
complex. Note that the van der Waals coefficient Cg ^ 
n^^ strongly increases with the atom's principal quan- 
tum number n. For typical values rt 60 this yields an 
enhancement of about 10^^ compared to interactions in 
low-lying states, implying a drastically different excita- 
tion dynamics. On the one hand, the resulting level shifts 
lead to a strong suppression of Rydberg excitation [l5|. 
due to an excitation blockade of close atoms [l^. On 
the other hand, the interactions perturb the atomic dark 
states 17j_18|, and admix the dissipative intermediate 
state [S, 11|. To account for these effects, we start from 
the von Neumann equation 
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for the TV-body density matrix p^^'' of the gas. The Lind- 
blad operator C describes spontaneous decay of the in- 
termediate state, but can also include finite laser band 
widths, denoted by 712 and 723, respectively. 



To obtain the steady state populations of eq.(l3]), we 
transform it to a many-body rate equation. For clar- 
ity we start from the simple case Cg = for which the 
A'^-particle density matrix factorizes, and eq.(I3]) reduces 
to a set of single-atom optical Bloch equations. Upon 
adiabatic elimination the dynamics of the atomic level 
populations follows 
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where the coefficients are straightforwardly obtained 
as a function of the laser parameters. 

This is a common approximation to the long-time dy- 
namics of two- level systems fl^ . Applications of the de- 
scribed elimination procedure to three-level atoms are, 
however, rather scarce, since it often leads to negative 
transition rates [2fl| . To resolve this obstacle we propose 
a linear transformation that removes the negativity and, 
at the same time, preserves the correct steady states of 
the underlying von Neumann equation ([3]). For fii < 1^2 
this is accomplished by adding a correction matrix 



(T2i(l - ^21) + (^32 

Aa = I -0-21 + 0-31 - (T32i?32 -0"2l(l — i?2l) - Cr31-R31 
(T21 — 0-31 + (T32i?32 — 1732 + 0-3ii?3i 
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0'2li?21 + 0-32(1 - RZ2) 
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to the original coefficient matrix a in eq. Q , where Uap — 
(lact/sl — 0'a(i)l^ and Rafi are free parameters. With this 
definition of the any negative rate coefficient aafi is 
set to zero, while the additional terms compensate for the 
according changes of the steady states. Consequently, the 
parameters Rafi are chosen such that the transformed 
rate equations p^'^ = (a + Aa)p'^*' yield steady states 
identical to those of the original equation p^'^ 
Explicitely, this condition gives 
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a22 (a3i + ass) + 032 (021 - a23) 

for the transformation coefficients in eq.®. 

Having established a proper rate equation description 
for the single particle dynamics, the approach can be ex- 
tended to interacting atoms. As shown in [21'], this is 
straightforwardly accomplished by replacing the upper 
detuning of each atom by = A2 — Uij . Here, 

the sum only runs over atoms in the Rydberg state |3), 
such that the local detunings and, hence, the individual 



atomic transition rates become dynamical variables that 
depend on the entire many-body state of the system. As 
a result one obtains a many-body rate equation where 
the dynamics of the ith atom is governed by its indi- 
vidual transition rates d^j^ which depend on the actual 
Rydberg atom configuration through the local detuning 
A2*''. Although this rate equation still covers the expo- 
nentially large number of all 3^ many-body states, it 
can be efficiently solved via classical Monte-Carlo sam- 
pling. Starting from the initial state with all atoms in 
their ground state \a) — |1), the steady state is obtained 
by performing repeated random transitions according to 

the probabilities p^'L.^ = 5t ■ a^*^ to make a transition 
within a time step St. In this way, we are able to ob- 
tain the fully correlated A^-body state-distribution of the 
particles, and in particular the average level populations 

Paa = N~'^ J2i Paa- SiuCC 
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the Monte-Carlo approach also allows to determine the 
imaginary part of the complex optical susceptibility 
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FIG. 2. Calculated probe beam absorption spectra for a 
Rubidium Rydberg-EIT medium at various densities. The 
lines show results of the Monte-Carlo simulations compared 
to quantum calculations based on a reduced-density matrix 
expansion 11] (symbols). The atoms are resonantly (A2 — 0) 
excited to 55S Rydberg states. The probe and coupling beams 
have Rabi frequencies of f^i = IMHz and ^2 ~ 2MHz, respec- 
tively, with linewidths of 712 = 723 ~ lOOkHz. 



where ^12 denotes the dipole moment of the probe tran- 
sition and eo is the permittivity of vacuum. In addition, 
the resonant real part of xi2 can be obtained from the 
nonlinear absorption spectrum using the Kramers-Kronig 
relations. 

In the following we consider the specific case of a Ru- 
bidium Rydberg gas. The atoms are excited by the 
probe laser (ili) from the |5S'i/2) = |1) ground state 
to the I5P3/2) = |2) intermediate state, while the cou- 
pling laser resonantly (A2 = 0) drives the transition 
between |5P) and a \nS 1/2) = |3) Rydberg state with 
ri2 = 2MIIz. We include the intermediate state decay of 
7 = 6.1MHz and further assume realistic linewidths [§] 
of 712 = 723 = lOOkHz for both beams. 

Fig [2] shows the calculated absorption spectrum at dif- 
ferent densities for ili = IMHz and n = 55. In order to 
check our results at low densities we have also performed 
simulations based on a reduced-density expansion of the 
von Neumann equation ([3]), as described in [ll|. The 
excellent agreement between these entirely different cal- 
culations attests to the quality of both approaches. At 
the lowest density, Rydberg-Rydberg atom interactions 
are ineffective, giving a small resonant absorption due 
the finite laser linewidths 712 and 723. At higher den- 
sities the interactions lead to a significant suppression 
of the resonant transmission. In accord with recent ex- 
periments Q, the position and width of the absorption 
minimum are, however, largely unaffected by the interac- 
tions, which are as large as C6(47rp/3)^ « 20MIIz at the 
highest density of 5 • lO^cm"'^. Apparently, this is due to 
the van der Waals blockade that prohibits simultaneous 
Rydberg excitation of close atoms. 

The density dependence of the complex susceptibility 
is shown in FiglJl For low and moderate densities our 



Monte Carlo calculations and the reduced-density matrix 
expansion 11] give consistent results for both the real and 
imaginary part of pi2 (and hence of Xi2)- As expected 
from a low-density expansion the latter tends to deviate 
with increasing densities, with significant deviations oc- 
curring only at rather high densities > lO^^cm^'^. As 
we show in the following, the Monte-Carlo calculations 
yield the expected high density limit, and should thus 
applicable for arbitrary densities. 

As the density increases the imaginary part of X12 
starts to saturate while the real part develops a maxi- 
mum and goes to zero at high densities. In fact, this 
behavior can be understood from simple arguments. At 
very high densities, a single Rydberg atom strongly shifts 
a large number of surrounding atoms out of resonance. 
As a consequence the Rydberg laser appears far detuned 
for the majority of atoms, such that they act as an effec- 
tive two-level medium. Taking the limit 00 the 
asymptotic high-density steady state value of pi2, hence, 
approaches 
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where 512 = (7 + 7i2)/7i 523 = (27 + 3723)/(7 + 723)- 
This simple limit is indicated by the horizontal arrow 
in Figl3] and gives very good agreement with the Monte 
Carlo simulations. 

It thus appears reasonable to employ the derived high- 
density limit as a natural scale for xi2- In addition we can 
use the interaction blockade of Rydberg excitation [lB| to 
rescale the gas density, i.e. the abscissa in FiglJ] Com- 
paring the Rydberg atom density pryd — PaN^^ J2i P33 
obtained from the Monte Carlo calculation to the corre- 
sponding value p^ylj for vanishing interactions one obtains 




FIG. 3. Complex susceptibility X12 as a function of the gas 
density po under EIT conditions, Ai = A2 — 0. The reduced 
density matrix calculations (symbols) start to deviate from 
the Monte Carlo results (lines) for densities > lO^^cm"^. Re- 
maining parameters are identical to those of Fig[2] 



4 



0.025 -(a) 
0.020 - 



■ ■■■■ 



0.015 - 



0.000 



«»,5,«tl,»,[»* I 



10' lo" lo' 10.'" lo" 
Po [cm ^] 




10 10 10 10 10" 10" 
/bl 



FIG. 4. (a) Nonlinear absorption coefficient Im[xi2] as a func- 
tion of the density of Rubidium atoms for 0,2 = 2MHz and 
7i2 = 723 = 0. The different symbols correspond to differ- 
ent principal quantum numbers and probe Rabi frequencies: 
0,1 = IMHz, n = 50 (squares), Oi = 0.5MHz, n = 50 (cir- 
cles), Oi = IMHz, 71 = 70 (triangles) and Oi = IMHz, n = 50 
(diamonds), (b) After proper scaling [cf. eas. (|10|l and ([9}] all 
data points follow the simple universal scaling relation eg. ljlip 
(line). 
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of suppressed Rydberg excitations. Re-expressing the to- 
tal atomic density through eq. pH]) and scahng the probe 
beam absorption by eq.([Hl), yields a universal dependence 

of Xi2 = Xi^lX]^^ on. the laser parameters, atomic den- 
sity and interaction strength. We have verified this be- 
havior for a wide range of parameters, some of which 
are exemplified in FiglUfor 712 = 723 = (FigH^.)- In- 



deed, all data points collapse on a single curve (Fig|4lD) 
described by 
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This simple formula nicely illustrates the effects of excita- 
tion blocking on the optical susceptibility of the Rydberg- 
EIT medium: For /bi < 1 one finds a nonlinear absorp- 
tion proportional to the blockade fraction, which, how- 
ever, saturates at the two- level limit Xi2^'' for /bi ?J 1, 
i.e. at the onset of a strong excitation blockade. Since 
/bl ~ Vl\ for weak Rydberg excitation 22] this implies a 
finite cubic nonlinearity, as observed in [y] . 

In conclusion, we have presented a numerical method 
that permits to explore the optical response of an KIT 
medium in the presence of strong atomic interactions. In 
particular, we considered, Rydberg-Rydberg atom inter- 
actions, which yield strong photon-photon interactions, 
manifested in a nonlinear f2i-dependence of the optical 
susceptibility. At small atomic densities our Monte Carlo 
results are in excellent agreement with reduced-density 
matrix calculations, but in addition are shown to be ap- 
plicable over a much wider range of densities and inter- 
action strengths. Beyond the present calculations, the 
described approach permits efficient simulations of very 
large atom numbers, which in future studies should en- 
able detailed comparisons with experiments. Besides illu- 
minating the effect of Rydberg interactions, the revealed 
universal behavior of the probe absorption may be of use 
in experiments, by permitting quick estimates of optical 
nonlinearities and allowing to determine the medium ab- 
sorption from pure population measurements. While we 
have primarily focussed on nonlinear absorption phenom- 
ena the developed approach can also be used to study the 
optical response of far-detuned FIT media and, thus, to 
explore prospects for realizing strong coherent photon- 
photon interactions. 
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